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Abstract 



Theory of strong decays defines in addition to decay widths, also 
QJ . the channel coupling and the mass shifts of the levels above the decay 

■ thresholds. In the standard decay models of the ^Po type the decay 

vertex is taken to be a phenomenological constant 7 and such a choice 
leads to large mass shifts of all meson levels due to real and virtual 
decays, the latter giving a divergent contribution. Here we show that 
taking the microscopic details of decay vertex into account, one obtains 
new string width coefficient, which strongly suppresses virtual decay 
contribution. In addition for a realistic space structure of the decay 
vertex of highly excited states, the decay matrix elements appear to 
be strongly different from those, where the constant 7 is used. From 
our analysis also follows that so-called flattening potential can imitate 
the effects of intermediate decay channels. 

1 Introduction 



There is large variety of single- channel models, proposed decades ago, which 
describe spectra of hadrons with reasonable accuracy |lj . The most popular 
and widely used is the relativistic quark model of N.Isgur and coworkers 
for mesons [2] and baryons [3j, where effective constants are used for quark 
masses (constituent masses), as well as an overall negative constant (C < 
0), and several additional parameters for spin-dependent interactions. For 
heavy quarkonia the Cornell model based on nonrelativistic Schroedinger 
equation and linear plus Coulomb potential, was extensively exploited. 
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Most of the models proposed are rather successful in predictions of low- 
lying hadron masses and the idea, that relativistic quark Hamiltonian with 
confining and the gluon-exchange potential, can be derived from QCD, seems 
to be realistic. It was indeed done in Ref . |5] , using the Wilson loop and field 
correlator technic, where for quarks at the ends of the rotating QCD string 
the relativistic string Hamiltonian (RSH) was derived. The RSH contains 
several improvements over old models: 

i) At small L (low rotation) it reduces to standard relativistic quark 
Hamiltonian [H 121 |5l |6] , but with current quark masses used instead of phe- 
nomenological constituent masses. The resulting hadron masses calculated 
are exprssed through the former and the string tension a [51 EJ [7] . 

ii) The overall negative constant is absent while for a given quark the 
universal negative self-energy correction appears, calculated via a |8]; its 
presence is crucially important to reproduce linear behavior of the Regge 
trajectories. 

iii) At high L due to string rotation term, which naturally appears in 
RSH [5j and is absent in quark Hamiltonian models [H El El S], the Regge 
trajectories with correct slope and intercept are calculated P fTU] 

As a result, one obtains the formalism, derived from QCD with minimal 
number of the first-principle parameters (current quark masses, ct^, and string 
tension a; connection of two latter was found in [TT]). 

Theoretical calculation of hadron masses with the use of RSH in single- 
channel approximation was successful for all states below open decay thresh- 
olds (see [T2j and ^13j for charmonium and bottomonium, [14j for heavy-light 
mesons, [9l [TOl [15] for light mesons, and |T6] for higher pionic states). How- 
ever, for states above threshold RSH gives somewhat higher masses and one 
can expect that taking coupling to decay channels into account one obtains 
mass shifts of these levels down, closer to experimental values. 

To this end the channel-coupling (CC) models were formulated in [1], 
[T71 [18] . They are based on the presumed form of the decay Hamiltonian, 
which is usually taken to be the ^Pq model [20] . More forms have been 
investigated in |2T], with a conclusion that the so-called sKs model yields 
results close to that within the ^Pq one. Influence of the CC effects on the 
spectrum are significant and can be divided in two parts. 

First, the effect of close-by channels, when the energy of the level in 
question is not far from the two-body threshold (e.g. ?/'(3770) in connection 
with DD,DD* thresholds). As was found in [I71[22|, the overall shift from 
the sum of the nearest thresholds (e.g. for charmonium) is of the order of 
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(100-200) MeV. 

Another part of the mass shift is associated with the contribution from 
higher intermediate thresholds (e.g. of a pair of higher D and D mesons) 
and in this case convergence of such terms appears to be questionable. This 
topic was investigated in [23] and in the first paper of p3] the authors have 
introduced additional form factor for quarks to make the ^Pq vertex nonlocal 
and ensure convergence of the sum of contributions over thresholds. 

Therefore the structure of the string-breaking vertex becomes a funda- 
mental issue and one should try to find its properties from the basic QCD 
Lagrangian, which takes into account both confinement and chiral symme- 
try breaking. Such the strong decay Hamiltonian was derived from the first 
principles in [21] (which also supported by the sKs model in its relativistic 
version), the interaction kernel being simply the confining potential between 
the newly born quark (antiquark q) and original (possibly heavy) antiquark 
Q (quark Q). This constitutes the strong decay term in action of the form 
(in the local limit cf [2T] ) 

Seff = J d^Xtpg{x)M{x)tpg{x), (1) 

A^(x,xq,Xq) = a(|xg-XQ| + |xg-XQ|). (2) 

From ([1]), (|2]) one obtains the decay matrix element between the origi- 
nal state {QQ)n and decay products - two mesons {Qq)n2 and {Qq)m with 
relative momentum p, 

(3) 

Here the factor 7/123 accommodates spin-angular variables and the func- 
tions in on]) refer to radial dependencies only, r = c(u — v), c ~ 1 for heavy 
QQ masses. In the way it was derived in |24j , the M. (x, u, v) refers to the 
string between positions u of quark Q and v of antiquark Q, which breaks at 
the point x somewhere between u and v. It is clear, that the point x should 
lie in the body of the string, i.e. within the string width d from the axis 
of the string, see Fig.l. This implies the necessity of an extra factor in ([T]), 
©string (x, I), which is proportioual to the energy density of the string with a 
fixed axis 1 (the vector u — v in our case). 
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Figure 1: String breaking (pair creation) point x and heavy-light radii R2 
and i?3, shown together with the string of radius d between charges Q and 



Now the string density was studied both analytically [25] and on the 

lattice [26l |27] . In the field correlator method [28l EH] the string width d is 
proportional to small vacuum correlation length, A ~ 0.1 fm [HI [29], and 
therefore it is also small, d < 0.3 fm, for not highly excited hadrons. 

The string field density was computed in [251 ISSl [S^ and one can visualize 
there the field distribution in the string, exponentially decreasing far from the 
string axis. In lattice calculations similar estimates hold, but they depend 
on the way of probing the string fields: in case of a connected probe one has 
c^con ~ 0.3 fm [26] and in the case of a disconnected probe d is smaller, ddisc ~ 
0.15 fm [27] . A simple look into the configuration of large closed Wilson loop 
for the string and a smaller one for qq closed trajectory, shows that c?disc is 
closer to the string breaking situation. In what follows we shall take d to be 
somewhere between the two (lattice) values. In next Section we shall study 
the effect of the decaying string width, called the factor 6string(x, I), on the 
decay matrix element and resulting mass shifts of energy levels. 

2 The width-of-the string correction in the 
string-breaking action 

In [21] it was shown that the effective action of the qq pair emission in the 
field of static charges QQ, placed at fixed points, can be written as 



Q. 
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where the mass operator Ai {x, y) is to be found from the nonhnear (integral) 
equation 

M{x, y) = (- i -A J{x, y) (5) 

and the kernel J{x,y) = Jq^x.y) + Jqi^x^y) accounts for the fields in the 
string. Taking into account only colorelectric fields of scalar confining corre- 
lator D{x), one can present Jq^x^y) as 

rx ry 

jQ{x,y) = —{Ai{x)A^{y)) = dui dviD{u-v). (6) 
Nc Jq Jq 

Here Q is the position of the closest static charge Q{ or Q) in 4d space 
and the analogous term appears for the anticharge Q (or Q). Note, that 
in [21] it was tacitly implied that in the averaging is over the vacuum 
configurations, and the points x, y can be anywhere in the space, surrounding 
static charge. It is the property of the kernel J(x, y) that it is asymptotically 
large for coUinear x||y , but the direction of this vector can be arbitrary. 
That was enough for the proof of Chiral Symmetry Breaking (CSB) due to 
confinement, but in our case one needs a further specification. 

Namely, at the moment of creation the created pair qq must lie on the 
minimal surface of the Wilson loop of static charges QQ, i.e. on (or inside) 
the string connecting static charges. This means that we must replace in 

{A4A4) by {A4{x)A4{y)) string, where the latter acquires the string profile 
factor 0string(a;, y), proportional to the string density of colorelectric fields, 

estnng(x,l/)=e(x,l)e(y,l). (7) 

Here 1 = XQ(t) — XQ(t) is the string axis vector. For long string, |1| ^ A, 
one expects that ^ depends only on the distance x^ from the string axis, e.g. 

2 _ I(x-xq) X Ip 

x± — p • yo) 

To simplify matter and for rough estimates one can take ^ as a Gaussian 
function of distance to the center of the string, so we take 

e(x,I)-exp(v(x-^^4^)'), (9) 

where P ~ ^ ~ 0(1 GeV). 
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Insertion of QstTing{x,y) in its local form, 9string(x, x), into Q is easily 
integrated and yields for intermediate mesons with almost equal radius, R2 ~ 
i?3 (corresponding SHO parameters (32 = /^s, (3i ~ l/Ri) 

Here 

One should note that the expression (ITT!) is valid as an asymptotic esti- 
mate for large distances R2, R2 ^ d. Besides the approximation (Q does 
not take into account an additional suppression in the case of short strings, 
|1| <R2 + Rs- 

Thus ri{(3,p) plays the role of the suppression factor for high excited 
intermediate mesons. Indeed, radii of high excited mesons are growing 
with radial and orbital quantum numbers {n, L), R2{n, L) ~ ^J~L,^fn, and 
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3 Study of the decay vertex 



The string decay vertex, derived in Ref. [21], has the form ([2]) in the local 
approximation. In the standard ^Pq approach [20] it is assumed that one can 
effectively replace the kernel A1(x, xq,xq) by a constant. The same type of 
approximation was used in [301 EB E2] and also in [21] , where results were 
compared with the analysis of decays of ^'(3730) — )■ DD and T(4S') BB in 
[33j . Below we shall study the reliability of this replacement and show that 
the replacement of the kernel by a constant is not always valid, especially 
for high excited states. 

To illustrate this statement we consider the decay matrix element ([3]) 
with kernel M, written in momentum space, where the wave functions of 
heavy-light mesons are replaced by gaussians \l>("2,n3)(^g~) _ (^^^^ e-9^/2/3| 
(/32 = 0.48 for D-mesons and (32 = 0.49 for S- mesons): 



'nn2n3\fj a/AT /3| J ^ ^^'^ QQ 



Jnninaip) = ^^"L"" / ?/l23^S(cp + q)e 
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Here $(a; b; z) is the confluent hypergeometric function, 



<^{a\ b: z) 



a a{a + 1) z"^ 
^+5^+ 6(6 + 1) 2! 



(13) 



the wave function \l/Q^(cp + q) can by expressed as a series of oscillator wave 
function (see Appendix 2 of [32j for details). 

Another expression, which should be compared with (JT2i) . can be obtained 
by replacement of the kernel Ai to some constant iW^: 



/^"^ (d) 



?/i23^S3t(cp + q)e-''^/^^^ 



(2^ 



|3' 



(14) 



We consider nS {n ranges from 1 to 5) charmonium states {"^qq in ( IT2|) 
and f lTi|) ). while the flnal states are (or DD* , D*D*) in all cases, then 
1/123 is proportional to g^. It is of interest to notice that while for the IS", 25" 
and 35* states one can reproduce Ai by constant rather well, on the contrary, 
for the 45" and 55" states such the replacement does not work (see left parts 
of Figs. 2,3,4,5,6). The constants iW^ appear to be different in these cases: 



0.65 GeV, 0.8 GeV and 1.1 GeV for the 1^,2^ and 3^ states. 





p, GeV 



Figure 2: Proflles of decay matrix elements, Eq.(3) (scalar parts), for lS{cc) 
into DD (left panel) and lS{bb) into BB (right panel) calculated with decay 
vertex of Eq. (2), expression ( IT2l) - solid line, and for the constant decay 
vertex, expression (HM - broken line. 

Surprisingly, that for the nS bottomonium states the constant decay ver- 
tex reproduces results with the kernel Ai very well even for 45* and 55* states. 
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however constants are different for different bottomonium states: M^^ varies 
from 0.65 GeV for IS" state to 1.3 GeV for 55" state (see riglit parts of Figs. 
2,3,4,5,6), what we see in the case of charmonium states too. 




p, GeV p, GeV 



Figure 3: The same as in Fig. 2, but for the 2S{cc) into DD (left panel) and 
2S{hh) into BB (right panel). 




Figure 4: The same as in Fig. 2, but for the ?)S{cc) into DD (left panel) and 
3S{bb) into BB (right panel). 



4 Analytic and phenomenological study of un- 
quenched spectra 

Our RSH was derived in Ref. [5j starting from the Wilson loop for the qq 
system and using Nambu-Goto action for the corresponding string. In the 
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Figure 5: The same as in Fig. 2, but for the AS{cc) into DD (left panel) and 
4:S{bb) into BB (right panel). 



derivation presence of additional quark loops was neglected (quenched ap- 
proximation), basing on the l/N^ argument and additional (phenomeno logi- 
cal) numerical suppression of the quark-loop effects. It is the purpose of the 
present Section to study these effects analytically and phenomenologically, 
and compare them with lattice results in the forthcoming Section. 

The generating functional of heavy charges QQ, after integrating over 
other quark-loops, has the form 



DAexp CaWqq{A) det{mg + D{A)), 



(15) 



where Ca is the standard gluonic action and Wqq{A) is the external (fixed) 
Wilson loop of heavy quarks. The det term can be written in the path 
integral form 



det(mg + D{A)) = exp 



trln(lf^{D^z)e-^WUA) 



(16) 



where [D^z) is the path integration, s is the proper time variable, K = 

i -fo (^^^T^) ^'^^ ^qq{A) is the Wilson loop of sea quarks, while tr implies 
summation over flavor indices and space-time coordinates. 

Expanding in the number of sea-quark loops, one has the first correction 
term [2l]: 

DAexpCA (^Itr 1 ^D'z)e-^WUA)WQQ{A)'^^ . (17) 



^lloop 
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P, GeV 



p, GeV 



Figure 6: The same as in Fig. 2, but for the 5S{cc) into DD (left panel) and 
5S{bb) into BB (right panel). 

Integrating in (fTTIl over DA, one obtains the effective one-loop partition 
function, 

^iioop = -- / -{D'z)e-^x{W,,, Wqq), (18) 

where x is connected average of two loops, 

X ^ {W,,iA)WQQiA)) - {W,,iA)){WQQiA)). (19) 

Properties of x were studied in [251 [26] , where it was shown that one can find 
a simple expression for x for small distances between minimal area surfaces 
of both Wilson loops, 

X = ^eM-^SA) (20) 

c 

and is the minimal area of the surface connecting contours Ci and C2 of 
Wilson loops Wqq and Wgg, respectively, as shown in Fig. 7 One can see in 
Fig. 7 that the width of the bands in 5'a along time direction is of the order 
of i?2,-R3, where R2,R3 are the radii of intermediate mesons (Qq) and (Qq). 

Several properties of the effective partition function (ITS]) can be derived 
immediately: 

1) The general property of the unquenching process: since the integral 
/ ^ is obviously diverging at small s, one needs a renormalization step, which 
means that the string tension a in ( l20l) is the renormalized (by unquenching) 
version of the quenched a. 
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Q R3 



Figure 7: Connected average of two Wilson loops, expressed via area law 
for the surface S'a between two loops. Radii R2 and of two intermediate 
mesons {Qq) and {Qq) are also shown. 



Expanding averaged static Wilson loop with sea quarks, one obtains 

{WQQ\,ii = {Wqq) + 1-^{{WqqW,,)) + ±{{WqqW,,W,,)) + (21) 

where the bar over Wgg implies averaging over all paths, i.e. all contours 
of light quarks Cqq with the weight defined by the Fock-Feynman-Schwinger 
path integral, 

Wq,^ / ^(L»S)e-^iy,,(C,,). (22) 

J So S 

Correspondingly, one can write each Wilson loop and their products as (omit- 
ting correction terms independent of T) . 

{Wqq) = eM-VQQ{R)T) (23) 

{WQQWgg)) = eM-V^{R)T); (24) 
{^QQ^^qWrn)) = exp(-yAA(i?)T). (25) 
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Therefore the resulting QQ interaction appears to be dependent on T. Since 
Vqq{R) = a,,^R + Vge{R), while V^AiR) < Va{R) < VQ^iR) for large 
enough R, with increasing T the QQ system will pass from the purely con- 
fining regime Vqq to one-loop regime Va and then to two-loop regime Vaa 
etc. In the next Section we shall show that this type of transition was indeed 
observed on the lattice. 

As to the form of Va(-R), Vaa(-R) etc., one expects that for R < R2 + R3, 
where R2,R3 are radii of lowest (Qq), (Qq) states, the form of Va(-R) does 
not change, i.e. 

Va{R) ~ areuR, i? < i?2 + i?3 ~ Ifm. (26) 

In case of the cc system R2 = R3 = Rd ~ Rd* ~ 0.6 fm. The same is 
true for bb system with Rb ~ 0.5 fm. 

In a similar way one can treat Vaa and higher loop terms. As a result 
one can predict that the static potential can be defined from the sum (^T^ in 
the T- independent way for i? < 1 fm, 

Vst,tUR)-VQQiR)^a,,,R, R<lim. (27) 

For i? > 1.2 fm the situation is complicated and static (T- independent) 
potential cannot be defined in the strict sense, as was discussed above. In 
this case another approach can be used, namely, the expansion of the con- 
nected averages {{WqQWqq)) in the series over intermediate heavy-light me- 
son states, as was done in [I8j, and it is equivalent to the expansions in [1], 
dZI, [201 ED E21 in this way instead of Va{R) one defines the energy- 
dependent nonlocal interaction 

V (n n' F^-S- [ "^'P ^»2n3(q-p)^,[,n3(q'-P) .^oN 
\/i2i(q,q,i^j - 2^ / -—7^ 7-7 , (28j 

where subscripts 1,2 refer to the channels QQ and {Qq){Qq), respectively, 
while n2,n3 are quantum numbers of the mesons (Qq) and (Qq) with the 
wave functions ipn2 ^-nd ipna, and 

Xn^n^ir) = ^^^^^^ / ^e^qV.,(q)^„3(q)- (29) 

In a similar way instead of Vaa(-R)5 one defines the interaction V131 due to 
three-meson intermediate states. 
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As a result, the total Hamiltonian has the form 



H = Hk^n + Vqq{R) + Vi2i{'R, R', E) + Vi3i{R, R', E) + 



(30) 



As one can see, in fl5U]) the T-dependent interaction of fl25]) . (^^, fl2^ is 
replaced by the energy-dependent nonlocal interaction. 

Let us underline general properties of the new Hamiltonian (for an earlier 
discussion see [M]). 

i) For energies below all thresholds the interaction J2n>2 ^ini(R, R, E) is 
negative, which implies attraction on average from all higher intermediate 
states. Hence the linear potential in Vqq{R) is modified (flattened) by in- 
clusion of intermediate states. This attraction also persists in some energy 
region above thresholds, where the real part of Vini is still negative. 

ii) Due to strong reduction of overlap integrals of the type | |\l/fc(R)Vi„i 
(R, R', i?)\l//(R') 1 1 (as was discussed in previous Section, it is due to the string 
width effect), the series J2n>2^ini is fast converging and therefore only few 
terms are important. 

Summarizing the effect of sea-quark loop on the QQ interaction, one can 
say that there is no energy-independent (or time- independent) universal local 
interaction which can describe the dynamics of QQ system in the unquenched 
case. If one tries to simulate the effect of quark loops on the static QQ 
potential, then it should be an approximate local interaction, which is close 
to linear potential (JrenR for R < 1 fm, and becomes softer (flattening) for 
larger i?, which can be approximated by making o"ren the energy - and R 
-dependent. 

Such kind of flattening potential was introduced in [IS] to describe high 
excitations of light mesons and used latter in [121 [131 [H] for higher charmo- 
nium states. 



Vqq{R) = a{R) ■ R- a{R) = ao 



1 - 7o 



exp(,/o^(r - Ri)) 



B + exp(y^(r - i?i)) 



(31) 



Here -Ri ~ 1.2 fm is the distance, where the string can decay into two mesons, 
aoRi+2MQ ^ 2MQq, = 0.19 GeV^. Putting 70 = 0.40, the best description 
of radial excited light mesons was obtained in [15] . The modified potential 
Vqq{R), taken from [15], is shown in Fig. 8. 

The resulting light meson masses, taken from [15], are compared with 
experimental data in Table 1. For heavy quarkonia the role of the flatten- 
ing is less important. Below in Table 2 one can see corresponding effect in 
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5 



linear 




r (GeV"') 

Figure 8: Modified potential with parameters given in Eq. (13T1) . For reference 
simple linear potential witli a = 0.19 GeV^ is also plotted. 
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Table 1: Light meson masses Mo(nL) (in GeV) from RSH with pure hn- 
ear potential and their masses Mo{nL) in potential with flattening. The 
corresponding flattening correction (5fiat (in MeV), the self-energy correction 
AsE = ^^^^ the string correction A^tr = ^^^^^^-jjr^-^ , and the resulting 

mass Mtot = Mq + Sse + ^str + ^ge (where 6ge is the correction from the 
gluon-exchange potential) are given in GeV. 



State 


7 = 


7 = 0.4 


5flat (MeV) 


As£ 


^str 


M = Mo 


nL 


Mo(nL) 


Mo(nL) 


flattening 








IS 


1.347 


1.335 


-12 


-0.510 





0.725 


2S 


2.009 


1.944 


-65 


-0.342 





1.503 


ID 


2.167 


2.122 


-45 


-0317 


-0.087 


1.662 


3S 


2.512 


2.300 


-212 


-0.274 





1.937 


2D 


2.615 


2.428 


-187 


-0.263 


-0.058 


2.052 


4S 


2.931 


2.569 


-362 


-0.235 





2.252 


3D 


3.006 


2.647 


-359 


-0.229 


-0.043 


2.322 


IP 


1.802 


1.777 


-25 


-0.382 


-0.074 


0.071 


2P 


2.328 


2.213 


-115 


-0.295 


-0.030 


0.068 


IF 


2.479 


2.402 


-77 


-0.277 


-0.113 


0.048 


3P 


2.766 


2.472 


-294 


-0.249 


-0.020 


0.065 


2F 


2.876 


2.606 


-270 


-0.239 


-0.083 


0.047 


4P 


3.146 


2.731 


-415 


-0.219 


-0.015 


0.062 


3F 


3.233 


2.821 


-412 


-0.213 


-0.064 


0.046 



charmonium levels, which is of the order of several tens of MeV for high 
excitations. 

Another types of flattening potentials were suggested in [35]. One should, 
however, be careful with the large R behavior of their flattening potentials, 
which is bounded from above, and therefore the quarks Q, Q can liberate 
themselves and this effect contradicts the physical picture in QCD, when an 
unstable hadron decays into hadrons, but not into quarks. 
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Table 2: Comparison of the single-channel and flattening potential results 
for S,P,D states of charmonium with existing experimental data. 



State 


SC 


flattening 


exp 


IS 
2S 
3S 
4S 
5S 


3.068 
3.678 
4.116 
4.482 
4.806 


3.066 
3670 
4.093 
4.424 
4.670 


3.067 
3.74(4) 
4.040(33^1) 
4.421(43^1) 
? 


IP 
2P 
3P 


3.488 
3.954 
4.338 


3.484 
3.940 
4.299 


3.525 (iPi) 
~ 3.93(5P2) 


ID 
2D 
3D 


3.79 
4.189 
4537 


3.78 
4.165 
4.475 


3.77(l3p)i) 
4.153(3) 



5 Comparison to other approaches 

Here we compare our string decay picture with lattice data and other ap- 
proaches. On the lattice the topic of string breaking and QQ interaction 
above inelastic threshold was actively explored during last decade (for the 
first attempts see [36] and [3lj). A way to determine the static potential 
in unquenched case was suggested in [38j and several spectra calculations, 
including sea-quark effects, were done in [52] • Recently careful studies of 
spectra of excited hadrons with open channels were published in [IQ] and 
|41j . where the importance of inelastic channels was stressed. The difficulty 
of existing lattice approaches is the lack of the proper definition of a reso- 
nance state, which actually belongs to the continuous spectrum and requires 
either the continuous density description or the use of the Weinberg Eigen- 
value Method, described recently in the last paper of Ref. [19]. The first 
approach is made possible by the use of the finite volume, when continuous 
states are discretized and the resonance is defined by the scattering phase 
EZ]- The second approach, to our knowledge, was never used on the 
lattice. As to precise definition of the resonance parameters on the lattice, 
from [ini HI] one can see that it needs a lot of efforts and is expected in the 
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Table 3: Comparison of calculated mass Mtot (in GeV) with experimental 
data. 



Mtot 


Theory 


Exp. 


Mcogith) 


M{p{lS)) 


0.749 


p(0.775) 


0.666 


M{p{2S)) 


1.519 


p(1.465) 


1.479 


M{p{3S)) 


1.937 


p(1900)? 


1.849 


M{p{4S)) 


2.252 


p(2150)? 


2.166 


M{IP) 


1.25 


ai(1230) 
/i(1282) 




M{2P) 


1.82 


ai(1647)?a2(1732)? 








/i(1815)? 
/2(2157) mixing for 3^P2 




M(3P) 


2.14 




M(4P) 


2.435 




1.65 


M{1D) 


1.66 


P3(1690),p(1720) 
mixing for 2^51 — l^Di 










M{2D) 


2.05 


P3(1990)? 


1.989 


M{3D) 


2.32 


P3(2250)?,p5(2330)? 
{3D - 2G mixing(?)) 


2.249 


M{IF) 


1.96 


04(2000), /4(2018) 




M{2F) 


2.24 


/4(2300) 




M{3F) 


2.50 







nearest future. 

It is worth saying that the QQ potential, calculated on the lattice, is not 
sensitive to the effects of the virtual sea quarks at least for distances -R < 1 
fm (for the latest calculation see |12]). This result is in agreement with our 
discussion of the structure of the unquenched Wilson loop in the previous 
Section. 
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